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Abstract 

Recently E-665 reported the first measurement of soft {En < lOMeV) 
neutron production in deep inelastic scattering (DIS) off nuclei. We report 
the first theoretical analysis of the data. We find that the observed cross 
section can be quantitatively explained as due to the final state interac- 
tions (FSI) of low energy nucleons {E^ ~ 200 — AOOMeV) produced in the 
elementary fiN interactions. We argue that the data indicate strong a sup- 
pression of the FSI's of fast partons (hadrons) in DIS at high energies, and 
that studies of the soft neutron production would provide a new sensitive 
probe of the dynamics of FSI's in DIS both at fixed target energies and at 
the HERA eA cohider. 
PACS numbers:25.30Mr 



The final state interactions (FSI) in deep-inelastic scattering (DIS) off nuclei are 
known to provide a unique way of probing the space-time evolution of strong interactions. 
However, observing these effects is challenging, especially if one wants to study their 
energy dependence. A number of experimental studies have focused on the production 
of leading hadrons in the current fragmentation region, see and references therein. 
The data indicate that at high energies the absorption of the leading hadrons is very 
small and hence most of the FSI occur at smaller rapidities. Indications of such FSI were 
reported by E665 in [|]. This poses a serious problem for the planned studies of FSI's at 
the HERA collider in its eA mode of operation (for review of the project see the report 
of the study group ||^) since the present detectors are not sensitive to hadrons produced 
in the proton (nucleus) fragmentation region. Hence we argued in that a study of soft 
neutrons may provide one of very few windows at HERA for studies of FSI's. Indeed 
such neutrons have energies ~ Ea/A {Ea is the energy of the nucleus) and very small 
transverse momenta. As a result the HI and ZEUS neutron detectors have nearly 100% 
acceptance for these neutrons 0. 

The production of soft neutrons with energies En < lOMeV in DIS off nuclei appears 
to be one of the very few observables which can be studied both at fixed target energies 
using standard detectors of low energy neutrons and at a eA collider using a forward neu- 
tron calorimeter [^0. The key question which we address here is whether soft neutrons 
are sensitive enough to the FSI's in DIS to be of use for such studies. 

The mechanism of the soft neutron production at intermediate energies in pA scatter- 
ing (kinetic energies, Tj„c < 0.5 — 1 GeV) is reasonably well understood - such neutrons 
are produced in the preevaporation and evaporation of the residual nuclear system left 
after interaction of the projectile with the target, see e.g. Refs. |Q. For a recent review 
see [^. In the case of heavy nuclei these neutrons provide the major channel of "cooling" 
of the residual system. 

In the first approximation the total neutron multiplicity is proportional to the number 
of nucleons knocked out from the nucleus in the initial (fast) stage of the interaction. 
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Recently the first measurement of tlie soft neutron yields in DIS was performed by 
the &665 collaboration |[TT|] at FERMILAB. Neutrons with energies En < 10 MeV were 



detected in high energy muon scattering off a number of nuclear targets: D, C, Ca, 
Pb. A relatively small average multiplicity of such neutrons, {Nn{A)) was observed. In 
particular, for the lead target where the data are most accurate 

{N^{Ph){Er. < 10 MeV)) = 5 ± 1. (1) 

As far as we know the only high energy data available in the high-energy data bases 
on the production of soft neutrons come from the ITEP (Moscow) experiment |[T^ for 
incident proton momenta of 1.4 < p^"^ < 8.5 GeV/c. We find that the shapes of the 
spectra for the overlapping energy range: 7.5 MeV < En < 10 MeV - are similar. How- 
ever, the neutron multiplicity is much higher for the proton projectiles (about a factor of 
5-6 larger for p — Ph collisions at Ep ~ IQGeV) which reflects a much larger number of 
wounded nucleons in pA interactions. In order to compare pA and ^A data, we calculated 
the number of neutrons produced per wounded nucleon ^^.^ \{aN) ^"J^p/e"^^'' ' since in 
Glauber type models without secondary interactions this quantity does not depend on 
the projectile |]T2[. This quantity is still about a factor of two larger for the proton case, 
indicating that secondary interactions are more important for the proton projectile. 



Average x for the trigger used in [11| is ~ 0.015. For such x nuclear shadowing for Ph 
is very small: atot = (^imp.appr - i^shad with - r] = (Jshad/ CTimp.appr. ~ 5%. Using the Gribov 
theory of nuclear shadowing which expresses the shadowing effect, (Jshad through cross 
section of coherent diffraction in 7*A^ scattering, and the Abramovskii-Gribov-Kancheli 
cutting rules one finds that nuclear shadowing leads to an increase of the cross section 
of the diffractive events by aghad] the decrease of the cross section of single interactions 
from aimp.appr. to (Jimp.appr. — ^o'shad^ and emergence of the cross section of simultaneous 



interactions with two nucleons equal to 2ashad , see |M for the recent discussion. Taking 



into account that the number of soft neutrons generated by the mechanism we discussed 
increases approximately by a factor of two for the events where two nucleons are knocked 
out, we find that the shadowing effects lead to a slight increase of average neutron 



multiplicity by a factor (1 — r])~^ ^ 1.05 for the kinematics of [^. due to the shadowing 
effects. Hence, we will neglect this small effect in the further analysis. 

It is generally assumed that in DIS high-energy hadrons are formed at distances 

''formation ~ ^2 V / 

from the interaction point (see e.g. 0]), in which ph is the hadron momentum and the 
scale /i^ is < 0.5 GeV^. Hence at momenta above 10 GeV they are formed beyond the 
nucleus and only hadrons with energies < few GeV are involved in FSI's with the residual 
nucleus. This picture is consistent with a very weak A-dependence of the leading particle 
spectra produced in DIS at high energies, see e.g. PJlI]. The most conservative assumption 
seems to be that only slow recoiling nucleons produced in the elementary lepton-nucleon 
DIS in the target fragmentation region reinteract with the nucleus. The spectrum of the 
recoil nucleons can be approximated at small x and not very large energies, where the 
triple Pomeron contribution is still small, as 

1 da^'+P^^+^ exp{-apl) 

oc (3) 



a^p dz(Ppt i/i 

with = (p2)i/2 _ o.4GeV/c, see discussion in [|12|. Eq.^ is consistent with the 

limited data of the BEBC neutrino experiment [|10|. Here z is Feynman x for the nucleons. 
Typical kinetic energies of the nucleons produced in the nucleus rest frame are ~ 200- 
400 MeV, and therefore by eq.2 they are formed within a modest nucleus radius. Soft 
neutrons should be produced both in the decay of the hole formed by the removal of a 
nucleon in the elementary lepton-nucleon DIS, and in the reinteractions of the secondary 
slow nucleons originating from the elementary DIS process. 

Contribution of these two effects should be considered as a lower limit for the rate 
of the soft neutron production. This is because there exist other processes for which 
the standard time formation arguments do not hold as well. For example the slow 
pion absorption process is known to increase the excitation energy of a nucleus by 60 
- 70 MeV and yield, on average, several additional soft neutrons per event, see e.g. P. 
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Contributions of of high-energy parton rescatterings, formation of hadrons inside the 
nuclei would further increase the neutron yield. 



We focus on the case of the DIS scattering off lead, since the muon data |11 



much more accurate in this case. To calculate the rate of the soft neutron production we 
considered the following model of the reaction: (i) a nucleon is removed from any point 
in the nucleus with a probability proportional to the nuclear density; (ii) an energy Wn is 
assigned to the produced nucleon based on cross section of the process fi+N fi'+N'+X 
as given by Eq.(3); the interactions of the nucleon produced in DIS are modeled with 
a Monte Carlo code for hadron -nucleus interaction which describes all stages of the 
process: the cascade, which includes knock-out of secondary nucleons, production and 
subsequent interaction of pions; preevaporation; and evaporation of neutrons and charged 
particles. This Monte Carlo code is very close to the codes are used over the years to 
describe cascades in pA scattering at En < IGeV, see e.g. HJ^. Note that in the original 
experimental paper |TT]] only a qualitative analysis of the data was presented, assuming 
that neutrons originate from the decay of the hole formed by removal of a nucleon in 
the DIS fiN interaction. Production of the soft neutrons in the FSI's of the nucleons 
produced in the primary DIS fiN scattering was not considered though in our calculation 
these FSI's provide most of the soft neutrons. 



The code was tested using available pA data at intermediate energies [|T5|-|T^. We 
found that the code produces a good description of the neutron yields in the kinematics 
of interest, see Fig. 1,2. The curves are about 20% above the data which is consistent 
with the typical 20% accuracy of such codes. This may also reflect uncertainties in the 
normalization of the data. 

We find that the calculated spectrum of the soft neutrons is consistent with the E-665 



data as reported in the erratum |Tl]], see solid line in Fig. 3. 

We want to emphasize that the calculated yield of neutrons depends weakly on the 
model used for the spectrum of nucleons produced in the elementary reaction. In par- 
ticular, assuming that all nucleons are produced with energy of 200 or 400 MeV (dashed 
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and dotted curves in Fig. 3) practically do not change our result. To illustrate further 
our weak sensitivity to the model of the nucleon production we present in Fig. 4 the mul- 
tiplicity of the produced neutrons for different cutoffs in En- One can see that it weakly 
depends on the kinetic energy of the generated nucleon, Wn, for the kinetic energies of 
interest: 200 MeV < Wn < 500 MeV. Note here that (Wn) for the model corresponding 
to Eq.| is ^ 300 MeV. 
We estimate 

{Nn{Pb){En < 10 MeV)),_ = 6 ± 1.5, (4) 

which is reasonably close to the experimental number of 5 ± 1. 

It is worth noting that there is a trend in the E665-data for {Nn{Pb){En < 10 MeV)) 
to fall with increasing qo = u (see Fig. 2 of |11[). This may be due to two causes: (i) 



the decrease of the probability for the leading hadrons to reinteract with nuclei with 
increasing u as observed in |jl|], and (ii) a softening of the spectrum of the nucleons 
produced in the elementary process with decreasing x (larger u in the data sample of 
E665 correspond to smaller average x). 

The multiplicity obtained in our "minimal" model of soft neutron production seems 
to leave very little room for the further processes of FSI of fast hadrons in the nuclei. The 
only alternative we could think of is that suppression of the FSI of produced nucleons 
starts at much lower energies than it is usually thought, say 1-2 GeV. 

To illustrate the sensitivity of the E-665 to a model of the FSI let us estimate {Nn{A)) 
in a class of the models where it is assumed that the leading quark can interact with 
effective cross section (Je// of 10-20 mb, see e.g. discussion in P,|l8|. The number of extra 
interactions due to this mechanism can be estimated as 

m = l^^eff I d'bT\b), (5) 
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where T{b) is the standard thickness function T{b) = J dzpA{b,z). In the case of lead 
target this leads to ~ 0.75cre// additional interactions per event where cXe// is measured 
in fm^. Low energy nucleons produced in these interactions would have the energy 



distribution close to the one in the elementary reaction and hence generate soft neutrons 
with the same efficiency as the "minimal" mechanism. Thus cTg// = 10(20)m6 would 
lead to the estimated soft neutron multiplicity of 8.75± 1.75 (12.5 ± 2.5) (to produce 
a conservative estimate we include here the pA data based adjustment factor of 0.8 for 
the theoretical value of eq.§ which is well above the experimental value of 5±1. We are 
planning to combine our MC code with a number of MC codes for production of high 
energy particles in eA scattering to constrain parameters of these models using E665 
data. 

It would be really challenging to reconcile a low neutron multiplicity observed by E665 
with the indications of a noticeable FSI's in production of faster particles in the nuclear 
fragmentation region obtained in the same experiment, see discussion in [n5|. This may 
require a more complicated dynamics of FSI's than currently envisioned. It also calls for 
a systematic experimental study of the energy dependence of the neutron yield starting 
from the incident energies available in the HERMES experiment and in the forthcoming 
COMPASS experiment. One can expect a gradual decrease of the rate of soft neutron 
production with increasing of u and Q^. It would be interesting to study correlations 
between the neutron multiplicity and the spectrum of leading hadrons (^-distribution, pt 
broadening, etc). It is necessary also to repeat the E665 experiment at higher energies 
to check rather amazing finding of this experiment of the low rate of production of soft 
neutrons. If the decrease with Einc and low neutron multiplicity at > 200 GeV are 
confirmed, the soft neutrons would provide a perfect tool to look for relatively rare final 
state interactions in DIS at the HERA collider at small x by selecting events with much 
larger than average neutron multiplicity. For example it would be possible to investigate 
the Pt broadening effects [pOjpIl via a study of the dependence of (A^„) on the transverse 
momentum of the leading hadron in the current fragmentation region. 

At the same time there are a few questions which could be addressed at the interme- 
diate energy lepton facility, like TJNAF. The first question concerns the estimate of the 
contribution to the soft neutron multiplicity from the decay of the hole state produced 
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by removing one nucleon in DIS. This problem is tightly connected with the study of the 
one-hole nuclear spectral function which is the subject of a number of (e, e'p) experiments 
that are under way at TJNAF. It would be of interest to measure the spectrum of neu- 
trons produced in the decay of the hole states, and its dependence on the missing energy. 
Hall C may provide possibilities for such study using a neutron detector in a combination 
with a high resolution (e, e'p) experiment. Such measurements also may be performed 
as a part of the recently approved (e, e'pn) experiment |Q. The experiment of the same 
kind but measuring yield of soft neutrons as a function of the energy transferred by elec- 
tron to nuclear proton (low energy resolution) can help to obtain the reliable estimate 
of the soft neutron multiplicity due to the FSI's of the intermediate energy protons with 
medium and heavy mass nuclei in lepton-nucleus interaction. 

We would like to thank K.Griffioen numerous suggestions and providing details of 
the E665 kinematics. We also thank L.Frankfurt, A.E.L.Dieperink and M.Sargsian for 
discussions. 
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FIGURES 




100 120 140 160 

E, MeV 

FIG. 1. Comparison of the results of the Monte Carlo cascade-evaporation calculation of 



the neutron spectra in p + Pb n + X process with the ITEP data |15| at Pp = lAGeV/c 
(lower curve) and Pp = 2GeV/c (upper curve) for 9n = 120°. The overall experimental errors 
which are not shown in the figure are ~ 20% 
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FIG. 2. Comparison of the results of the Monte Carlo cascade-evaporation calculation of 
the neutron spectra in p + Ph n + X process at £^"'^=113 MeV with the data |T^ . 
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FIG. 3. Comparison of the results of the Monte Carlo cascade-evaporation calculation of 



the neutron spectra in + P6 — > n + X process with E-665 data 
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FIG. 4. Multiplicity of neutrons produced in the process where a nuclcon with energy W 
was produced in Pb for different energy intervals of neutron energy. Dotted, solid, dashed 
curves are multiplicities of evaporated nucleons with kinetic energies T„ for < r„ < 6 MeV, 
< T„ < 10 MeV, < T„ < 50 MeV; dashed-dotted curve is the total nucleon multiplicity. 
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